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Abstract. The light curves of the two images of the dou- 
ble quasar Q0957+561 as obtained by Kundic et al. (1997) 
are almost identical, except for an overall time delay and 
scaling factor. This allows us to put limits on the amount 
of microlensing that took place during the time interval 
corresponding to the monitoring observations. We per- 
form numerical simulations in which we model the mi- 
crolensing behaviour of the (halo of the) lensing galaxy 
in the system. We test "MACHO-masscs" ranging from 
1CT 8 to 10 _1 M and quasar sizes from 10 14 to 3 x 10 15 
cm. Statistically comparing the expected microlensmg- 
induccd changes from 100000 simulated light curves over 
a period of 160 days with the (lack of) observed fluctu- 
ations, we can constrain regions in the parameter space 
of MACHO mass and quasar size with various degrees of 
confidence. In particular, a halo consisting of objects at 
the low end of our mass scale can be ruled out with high 
confidence for a small quasar size. A halo consisting of ob- 
jects with 10~ 2 or 10 _1 M Q cannot be ruled out yet, but 
it should produce MACHO induced fluctuations in future 
observations. We also test halos with only 50% or 25% of 
the mass in compact objects; constraints here are a bit 
less stringent. 
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1. Introduction 

It has been known for almost three decades that rota- 
tion curves of galaxies remain flat even beyond the visible 



matter (Rubin & Ford 1970 , Rubin et al. 1985 ). As an ex- 
planation, Ostriker et al. (1974) and Einasto et al. (1974) 
suggested that galaxies are surrounded by large halos con- 
sisting of "dark matter" and dominating the total mass of 
the galaxies. Over the years a very large number of possi- 
ble "candidates" for this dark matte r hav e been sug geste d 
(for recent reviews see, e.g., Bahcall 1997 or Raffelt 1997 ). 
They can broadly be divided into "elementary particle" 



candidates (e.g. massive neutrinos or axions) and "as- 
trophysical" candidates (e.g. black holes, brown dwarfs, 
" Jupiters" , comets) . Not much progress has been made to 
date in identifying the elusive dark matter despite major 
efforts in many directions. In fact, a new branch of physics 
established itself - "astro-particle physics" - whose major 
goal is the solution of the dark matter problem. 



More than 10 years ago Paczyhski (1986) proposed a 
direct test to prove or reject the possibility that the halo 
of the Milky Way consists of dark compact objects by way 
of gravitational microlensing: If the brightness of at least a 
million stars in the Large Magellanic Cloud (LMC) could 
be regularly observed, at any given time about one of them 
should be magnified due to a halo object (in the following 
referred to as "MACHO" , for "massive compact halo ob- 
ject") passing in front of it and focussing the light rays to 
the observer. The light curve of the affected background 
star should show a very characteristic and achromatic be- 
haviour. 

Soon thereafter various groups started big observa- 
tional programs to investigate this promising possibility, 
and the first microlensing events were found in 1993 (Al- 
cock et al. 1993; Aubourg et al. 1993). The most recent 



results indicate that about 501^0 % °f the dark matter in 
the Milky Way halo can consist of such objects (Alcock 

+0.3 
-0.2 



et al. 1997). Their most likely mass range (M « 0.5 



M Q ), however, is higher than originally suggested for gen- 
uine "brown dwarfs", but the uncertainty is large. The 
largest problem in these experiments is the small proba- 
bility for a microlensing event: the "optical depth" (the 
fraction of background stars that are significantly affected 
by microlensing) is less than 10~ 6 . 

In a very different optical depth regime, gravitational 
lensing can also be used to test whether halos of other 
galaxies are made of compact objects (Gott 1981 ): In cases 
of good projected alignment along the line of sight be- 
tween a background quasar and a foreground galaxy, the 
galaxy can produce multiple images of the quasar^. In this 



1 Note that here the distances involved are about five orders 
of magnitude larger than in the microlensing searches towards 
Send offprint requests to: R. Schmidt, email: rschmidt@aip.de the LMC (few Gigaparsec rather than some 55 kpc), and the 
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case the light bundles from the quasar pass through (the 
dark halo of) the foreground galaxy on their way to the 
observer. These light bundles probe the graininess of the 
halo: if the dark matter there consists of some kind of 
elementary particles ("smoothly distributed matter" on 
astronomical length scales), the light bundles should be 
unaffected; if it is made of compact astrophysical objects, 
the light bundles can probe it. The measured brightnesses 
of the quasar images should vary as a function of time due 
to the changing relative posi tions between lens, source and 
observer (Chang & Refsdal 1979 ). The surface mass den- 
sity (optical depth) in these cases is high enough (of order 
unity) to basically affect (change) the measured flux of a 
quasar image all the time. 

Since quasars are intrinsically variable as well, it is not 
trivial to decide whether an observed variability is intrin- 
sic or microlensing-induced. In the case of multiple images 
of one quasar, however, the intrinsic fluctuations should 
show up coherently in all the images - modulo a certain 
time delay - whereas the microlensing changes occur un- 
corrected in the various images. Once well-sampled light 
curves of two (or more) quasar images are obtained and 
the time delay At (due to the different light paths of differ- 
ent light bundles) and the magnitude difference Am (due 
to different magnifications) is known (both of which are 
not trivial), one can shift the two light curves in time and 
magnitude by the appropriate amounts At and Am and 
subtract them from each other. All remaining fluctuations 
in the "difference light curve" must be due to microlens- 
ing. In particular, if the difference light curve is flat, this 
indicates that there was no microlensing going on during 
the period of observation)^. 

The gravitationally lensed double quasar Q0957+561 
(Walsh ct al. 1979] ) has been observed/monitored for al- 
most two decades by many groups (see e.g., Schild & 
Thomson 1995) with the goal to measure the time de- 
lay and subsequently determine the Hubble constant. The 
time delay is now established firmly at At = 417 ± 3 



days (Schild & Thomson |1997j K undic et al. |1997[ sub- 
sequently K97; Oscoz et al. 1997). In addition, a number 
of gravitational lens models have been published for the 
Q0957+561 system which reproduce the observed (optical 



and radio) image configurations well (Falco et al. 1991 



Grogin & Narayan 1996), so that the parameters relevant 
for microlensing at the position of the quasar images can 
be determined from these models with an accuracy of a 
few percent. 

With both the time delay established and good models 
being available, the double quasar Q0957+561 can hence 
be used as a test for massive halo objects by comparison 



optical depths involved are of order unity (i.e. six orders of 
magnitude higher). 

2 The effect of "flat" microlensing - periods in which mi- 
crolensing produces a constant (de-) magnification - are taken 
care of by a different magnification ratio. 



of the light curves of images A and B. We do this here 
by simulating the microlensing effect numerically for MA- 
CHOs of different masses. We analyse the resulting mi- 
crolensed light curves and show how often microlensing- 
induced changes of certain amplitudes are to be expected 
for certain MACHO masses. Finally, we compare our nu- 
merical simulations with the well-sampled data set of K97 
and thus constrain the MACHO masses in (the halo of) 
the lensing galaxy. 

In Sect. 2 we describe briefly the monitoring data set of 
K97 for the double quasar and the simulations we perform 
in order to compare the observations with microlensing 
light curves due to different masses. In Sect. 3 we show our 
results for a MACHO mass range from 10~ 7 < M/Mq < 
10 -1 , and in the final Sect. 4 we present our conclusions. 



2. Observations and simulations 

We first briefly present the macro model of the lens system 
and the recent data set on the quasar Q0957+561 which 
we use. We then define the relevant numbers and param- 
eters and illustrate how the microlensing light curve de- 
pends on the MACHO masses and the quasar size. Finally, 
the numerical technique is described which we employ to 
produce the microlensing light curves and to analyse them. 



2.1. Lens model and optical data 

The double quasar Q0957+561 (z = 1.41) is gravitation- 
ally lensed by a galaxy (z = 0.36) and its associated galaxy 
cluster. The complex structure of the quasar in the ra- 
dio regime places strong constraints on theoretical models 



(e.g. Gorenstein et al. 1983, 1988), so that well determined 
values for the surface mass density and the local shear at 
the positions of the quasar images can be obtained (Falco 
et al. 1991, priv. comm., Grogin & Narayan 1996). 

The normalized surface mass density k = E/E cr it is 
the projected mass density of the lensing galaxy plus clus- 
ter along the line of sight, normalized by the critical sur- 
face mass density (Schneider et al. 1992) 



A, 



AttG D d D d 



= 0.92/ig/cm 2 



(1) 



Here D d , D s and D ds are the angular diameter distances 
between observer, deflector and source, respectively (the 
velocity of light is denoted by c and the gravitational 
constant by G). The underlying cosmological model is 
an Einstein de-Sitter universe with a Hubble constant of 
H = lOO/i.kms^Mpc- 1 . 

The local shear 7 represents the tidal field at the image 
position due to the matter outside the beam. We adopt 
values for the surface density of ka = 0.32 for image A and 
kb = 1-17 for image B, for the local shear at the image 
positions we used ja = 0.18 and 7^ = 0.83, respectively 
(Falco priv. comm.). 
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Recently, K97 presented the results of two years of well 
sampled monitoring observations of Q0957+561 with the 
Apache Point Observatory. With their g band data, they 
confirmed the value of 417 ± 3 days for the time delay At 
between the two lensed light paths and a best value for 
the magnitude offset between the (time-corrected) quasar 
fluxes of 

Am AB := < m A (t) - m B (t + At) > =0.118 mag. (2) 

In order to quantitatively estimate the effects or lim- 
its of microlensing on the light curves of Q0957+561 A 
and/or B, we determined the difference light curve be- 
tween the two quasar images in the following way. First, 
the flux-corrected and time shifted light curve of image 
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B (i.e.: m B (t) 



,(t) = m B (t + At) + 0.118 mag) was 



subtracted from the light curve of image A: 
Am(t) = mA(t) — m B (t), 



(3) 
light 



which measures the deviation between the two 
curves. 

Since the light curve had been sampled only for dis- 
crete points in time, we had to interpolate between the 
two closest points of light curve B before and after the 
instant of time in which a point of light curve A had been 
determined (or vice versa): 



Fig. 1. Difference light curve of images A and B (shifted in 
time and magnitude) of the quasar Q0957+561. To guide 
the eye, dashed lines are drawn at differences of +0.025 
mag and -0.025 mag. 

difference light curve and a horizontal line (i.e. the hypoth- 
esis: no detectable microlensing-induced changes within 
the measurement uncertainties) defined by: 



1 



N 



N 



j(tj) = Ci m' B (U) + Cfc m' B (t k ), 



(4) 



Am 2 {tj 



(6) 



where & 



t 



(tj - U)/(t k - fj) and ^ + Cfc = 1 and U < 



N is the number of data points of the difference light curve. 
This x 2_va l ue measures the goodness of representing the 



< ifc.; tj indicates the time at which data was taken for d a t a with Am = We obtain x 2 = 1 or 



image A, and U,tk are the closest times before and after 
ti for which data exists for image B. The "difference light 
curve" is defined as: 



Am(tj) = rriA(tj) — m B (tj). 



(5) 



The measurement uncertainties cr(tj) were added quadrat- 
ically for each combined pair of data during interpola- 
tion and subtraction. The resulting difference light curve 
is shown in Fig. |[ We also determined the difference light 
curve by interpolating the A-light curve (rather than B) 
and obtained a very similar result. 

We used 36 observations of image A and 40 observa- 
tions of image B out of the 99 available observations from 
the years 1995 and 1996 (K97) since we had to restrict 
ourselves to the about 160 days of "overlap" between the 
image A light curve and the time-shifted image B light 
curve. The time axis in Fig. [l] corresponds to the observ- 
ing epoch of image A. 

We detect no variation in the difference light curve 
with an amplitude greater than sa 0.05 mag. And consid- 
ering the error bars we find that the difference light curve 
would even be consistent with Am = mag. There is 
also no systematic gradient apparent in the data, which 
would be the signature of a long term microlensing event, 
produced by a relatively massive MACHO. The small vari- 
ance can be quantified. We determine a x 2 between the 



X 2 = 1.2, 

depending on which of the two light curve sets was inter- 
polated (see above). This means there is no statistically 
significant deviation between the light curves of the two 
images. 

However, in the difference light curve in Fig. |l| there 
is a small peak seen around day 750 (K97). We tried in 
many ways to establish the significance of these half dozen 
or so points which are slightly but coherently above the 
"zero" -line. Interestingly, if one goes back to the original 
individual light curves by K97, it can be seen that at the 
Julian date corresponding to the peak in light curve of im- 
age A, a similar peak was seen in the light curve of image 
B. This could indicate a systematic effect during the ob- 
servation. Nevertheless, this "peak" in the difference light 



curve seems to have been seen by Schild (1996b) as well 
with independent data taken at a different observatory. 

If real, such a peak in the difference light curve could 
be interpreted as a short-duration, small-amplitude mi- 
crolensing event in image A. In that case it could provide 
valuable information on the mass of possible microlensing 
objects in the halo of the lensing galaxy. Since we cannot 
prove with our data that this peak is real (rather than a 
statistical fluke), in this paper we will argue more conser- 
vatively that we do not detect any microlensing changes 
Am max — Am m ; n higher than 0.05 mag; we will use this ar- 
gument to put limits on the masses of the MACHOs in the 
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halo of the lensing galaxy. This leaves open both the pos- 
sibility that this possible event is a fluke or that it is real. 
If the latter turns out to be true, one can even draw some 
stronger quantitative conclusions on the masses of possi- 
ble MACHOs, rather than just exclude certain regions in 
parameter space. 

2.2. Simulating microlensing light curves 

The exact shape of a microlensed quasar light curve de- 
pends on 

— the direction of the (projected) relative velocity vector 
between lens and source 

— the (projected) relative positions of the MACHOs 

— the masses of the MACHOs and 

— the size of the continuum emitting region of the quasar. 

We do not and cannot know the exact positions of 
the MACHOs. Hence we will not be able to predict or 
explain an individual microlensed quasar light curve (in 
contrast to the "low optical depth" regime of microlensing 
of stars in the Milky Way or Magellanic Clouds). However, 
we can determine and analyse microlensed light curves in 
a statistical sense; in particular, we investigate here the 
distribution of total magnification variations. 

We first determine the two-dimensional magnification 
variations due to microlensing at different positions in the 
source plane with the ray-shooting technique; we follow 
light rays backwards through an arrangement of MACHOs 
randomly distributed in the plane of the lensing galaxy, 
with surface mass density and shear as given by the lensing 
model by Falco et al. (1991, priv. comm.) (see Sect. 



The density of the deflected light rays in the quasar plane 
corresponds to the relative magnification as a function of 



position (Wambsganss 1990). In these simulations, we al- 
ways use MACHOs with identical masses.^] 

We follow approximately of 10 10 light rays and col- 
lect them in the source plane in an array of 2500 by 2500 
pixels. We simulate microlensing light curves in the stan- 
dard way (see, e.g., Kayser 1986 or Wambsganss |l990| ) by 
evaluating the magnification along linear tracks across the 
magnification patterns with a physical length L, equal to 
the length the quasar traverses in the 160 day time span 
of the observed difference light curve shown in Fig. [I]. By 
assuming a projected quasar velocity relative to the mag- 
nification pattern of Vt = 600 kms" 1 (calculated using the 



method by Witt & Mao 1994 ), this length is given by 



L = 600 km s -1 x 160 days = 8.3 x 10 14 cm. 



(7) 



We have neglected any motion of t he MA CHOs relative to 
each other (Kundic & Wambsganss 1993 and Wambsganss 



In simulations with steep mass functions, e.g. Salpeter-like, 
most of the objects are near the lower cut-off; so the results of 
microlensing simulations with such mass functions are similar 
to th ose wi th all objects identical to the mean mass (Lewis & 
Irwin |l996|). 



& Kundic 1995). Since velocities of stars in galaxies are in 
general smaller than galaxy velocities, this effect cannot 
dominate the bulk velocities. It merely slightly increases 
the value of the transverse velocity of the quasar. 

2.3. MACHO mass and quasar size 

We calculated magnification patterns for both quasar im- 
ages with fractions of the halo mass contained in compact 
objects of 100%, 50% and 25%. For each of these mass 
fractions we produced magnification patterns of varying 
physical side lengths of 20, 200 and 2000 Einstein radii. 
The Einstein radius in the source plane is defined as (e.g. 
Schneider et al. |l992[) 



/4GM D S D, 



,1s 



1/2 



3.7 x 10 1 



l^ h -i/2, 8) 

Mq 



From this expression one can see that a single magnifica- 
tion pattern can be used to simulate light curves for var- 
ious MACHO masses because the physical length scales 
with the square root of the mass of the MACHOs§ the 
fixed length the quasar traverses in the source plane trans- 
lates into different numbers of pixels for different MACHO 
masses and pattern sizes. The only limitation we have for 
the investigation of various microlensing masses is the dy- 
namic range of the magnification pattern, which in our 
case is an array of 2500 by 2500 pixels. We used track 
lengths ranging from 10 pixels to 300 pixels, so that we 
could simulate three decades of MACHO masses with one 
pattern. With the three different side lengths, however, 
we were able to investigate the effects of MACHO masses 
M ranging from lO~ 8 Af0 up to 1O _1 M0 (in steps of fac- 
tor 10). Due to the "overlap" we could check some masses 
on magnification patterns with different side lengths and 
hence cross-check the results. 

It follows from these considerations that the quasar 
traverses more (fewer) characteristic lengths of the magni- 
fication pattern for smaller (larger) MACHO masses dur- 
ing the observation period of 160 days. This implies that 
the light curve is more (less) variable for smaller (larger) 
MACHO masses, so that one can derive limits on the MA- 
CHO masses from the microlensing variability. This is il- 
lustrated for three different mass scales in Fig. ||; the vari- 
ability of the microlensing light curves increases strongly 



4 An easy way to understand this "scaling argument" with- 
out invoking the concept of an "Einstein radius" is a follows: 
Consider a certain region in the source plane. Suppose that 
we know the surface mass density of the MACHOs in the lens 
plane. Let us raise the mass of all MACHOs by a factor q. 
If the surface mass density is kept constant in the lens plane, 
this corresponds to blowing up the length scale of the distri- 
bution of MACHOs by a factor ^/q. Since the deflection angle 
is proportional to the mass and inversely proportional to the 
distance to the lens, the magnification pattern in the source 
plane is thus also blown up by the same factor of ^/q. 
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Fig. 2. Magnification patterns and light curves for different MACHO masses. The three microlensing light curves on the 
right arc simulated by evaluating the magnifications along the white tracks in the three magnification patterns on the 
left. The track length is chosen in such a way that it corresponds to a 160 day light curve of Q0957+561 for an assumed 
transverse velocity of the quasar across the magnification pattern of i> t = 600 kms -1 . The magnification patterns are 
greyscale-coded from white (high magnification) to black (demagnification) and have side lengths of 4 Einstein radii 
(top), 40 Einstein radii (middle) and 400 Einstein radii (bottom), so that the (fixed) track length corresponds to 
MACHO masses of 10~ 3 M Q (top), 10 _5 M Q (middle) and 10~ 7 M Q (bottom). The microlensing parameters of the 
magnifications pattern are those obtained by Falco et al. ( |1991 , priv. comm.) for image A of Q0957+561; the surface 
density (in units of £ cr it) is 0.18, and the shear is 0.32 (aligned horizontally in this figure). The quasar is simulated 
with a Gaussian surface brightness profile with a half-width of 3 x 10 14 cm. 
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with decreasing MACHO mass. Qualitatively one can take 
from the magnification patterns in Fig. |^ that the proba- 
bility of observing no microlensing variation is practically 
zero for 10 _7 Mq MACHOs (bottom panel) whereas for 
the 10~ 3 M o MACHOs (top panel) this is not unlikely. 
These qualitative statements will be quantified in the next 
section. 

The size of the optical continuum region of the quasar, 
which we call quasar size in the following, has an impact 
on the shape of a light curve. For extended objects the 
total magnification can be calculated as a weighted mean 
of the magnifications at many points in the source plane. 
In practice the effect of the source size can be accounted 
for by convolving the two-dimensional magnification pat- 
tern with an appropriate source profile. Sharp features in 
the magnification pattern - especially the line-like caus- 
tics - are thus smoothed out by the brightness profile of 
the quasar. The amplitude (smaller for large sources) and 
the duration (longer for large sources) of variations in a 
microlensing light curve hence depend on the source size 
(see Wambsganss & Paczynski 1991 ). 

We adopted a Gaussian profile for the surface bright- 
ness profile of the quasar, where the source size is defined 
by the Gaussian width ctq. We used quasar source sizes 
ranging from ctq = 10 14 cm up to ctq = 3 x 10 15 cm (in 



steps of vlO). These quasar sizes are smaller than or of 
the order of one light day, which is an upper limit on the 
quasar size that was ob tained by Wambsganss et al. ( 1990 ) 
and also Witt & Mao ( |1994| ) for the quasar Q2237+0305. 
The case for even larger quasar sizes and smaller MACHO 
masses was e xami ned by Refsdal & Stabcll ( |l99l| , |l993D 
and Haugan ( 1996| ). 

We note that all our results for the MACHO masses 
and quasar sizes may be scaled with v = wt/600 kms -1 
for various values Vt of the transverse source velocity. MA- 
CHO masses scale quadratically with v because the scale 
of the magnification pattern is proportional to the square 
root of the MACHO masses. Similarly, the source sizes 
scale linearly with v. 

3. Results 



For each set of the two parameters MACHO mass and 
quasar size and for each of the two quasar images A and 
B, we analysed 100 000 randomly chosen tracks across the 
magnification patterns. The same tracks were used for dif- 
ferent source sizes. For each light curve we determined the 
difference between the highest and the lowest point of the 
light curve Am max — Am m i„. We call this quantity the 
total magnitude variation of the light curve. 

With these light curves, we can calculate the proba- 
bility p >c i of observing a total variation greater than or 
equal to some value d for each analysed parameter pair of 
MACHO mass and quasar size. As an example, in Fig. [| 
two integrated probability distributions p >( i are shown for 
MACHOs of mass lO^Mg and 10~ 5 A/ o . In these plots 



it is assumed that the quasar has a size of 10 14 cm and 
that the halo mass is completely made up of MACHOs. 
Three lines are shown per plot; the distributions for images 
A and B alone, as well as the joint probability distribu- 
tion where at least one quasar image has a total variation 
greater than d. One can see that much stronger microlens- 
ing variations are expected for small MACHO masses than 
for large masses on these short time scales. 

In the following, p >c [ refers only to the joint probabil- 
ity where the microlensing variation is observed in at least 
one quasar image. In order to compare the simulations 
with the observations, we calculated the joint probability 
P>o.05 for a grid of points in our parameter space of MA- 
CHO masses and quasar sizes. We inferred in Sect. 2.1 
(Fig. |l|) that the observations indeed do not show a to- 
tal variation greater than 0.05 mag, so that the quantity 
P>o.05 can be viewed as the confidence level at which we 
can exclude a particular parameter pair as not consistent 
with the observations. 

Various calculated values for p>o.o5 & r e given in Ta- 
bles [j] (for an assumed halo fraction of the MACHOs of 
100%), | (for a halo fraction of 50%), and | (for a halo 
fraction of 25%). The table entries for parameter pairs 
that are ruled out at a confidence level of 95% and above 
are highlighted in grey. The uncertainties given in brack- 
ets are the standard deviations from the mean of three 
different values for p>o.o5 that we derived from three inde- 
pendent realizations of each magnification pattern (using 
different random MACHO fields). The parts without en- 
tries are regions in the parameter space which could not be 
accessed because they were beyond the dynamical range 
of our simulations. 

We have illustrated the results from Table [j] in Fig. |[ 
In this plot the confidence levels are represented by the 
height and the colour of the plotted bars. The numbers 
from this plot and Tables 0, | and | show that MACHO 
masses in the region from 10 M down to 10~ 5 M can 
be ruled out to make up a sizable fraction of the halo 
mass in the lensing galaxy of Q0957+561 for quasar sizes 
smaller than 3 x 10 14 cm. In fact, for quasar sizes smaller 
than 10 14 cm, we can exclude MACHO masses down to 
10~ 7 M o . 

When looking at the uncertainties of the given proba- 
bilities in Tables [l| || and || we find relatively good agree- 
ment between the results from the three independent re- 
alizations. Only for the light curves with pixel lengths of 
300 pixels or 12% of the side length of the magnifica- 
tion patterns (corresponding to the smallest masses sim- 
ulated with each magnification pattern) the discrepancies 
between the probabilities get larger. The reason is that the 
size of the light curve becomes comparable to the whole 
field size and hence the light curves are not strictly inde- 
pendent of each other. In order to get statistically reliable 
results, one should use the p>o.o5 that was derived from 
the largest magnification pattern (in terms of side lengths 
in units of Einstein radii) that can be used to generate 
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Fig. 3. The probability p >( i of observing a microlensing variation greater than d in the quasar Q0957+561 either in 
image A (dotted line), image B (dashed line) or in at least one of the two images (solid line). In the plot on the left 
the MACHO mass is 0.1M©, whereas in the plot on the right it is 10 _5 M Q . In both plots the quasar size is 10 14 cm 
and it is assumed that the halo is completely made up of MACHOs. 




Fig. 4. "Exclusion" probability for certain MACHO masses. Three-dimensional visualization of the probabilities p>o.05 
(in percent) for measuring a total microlensing variation greater than 0.05 mag in a 160 day difference light curve of 
Q0957+561 for a particular parameter pair of MACHO mass and quasar size. The probabilities are indicated by the 
grey-shade of the bars (see the key), the relative scale is visualized by the bar height. The parameters of the blank 
field were beyond the dynamical range of our simulations. These probabilities are those from Table [I] for the largest 
available magnification pattern for each parameter pair. It is assumed that MACHOs make up 100% of the halo mass. 



the light curve sample in question. Incidentally, we find 
this way that some of the caustic patterns show coherent 
structures on scales as large as 200 Einstein radii. 



4. Discussion and conclusion 

In two years of observational data on the gravitation- 
ally le nsed double quasar Q0957+561 by Kundic et al. 
( 1997 ) no microlensing variation of the quasar larger than 
0.05 magnitudes was observed. In our microlensing sim- 
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Table 1. Probabilities p>o.o5 ( m percent) for measuring 
a total microlensing variation greater than 0.05 mag in a 
160 day difference light curve of Q0957+561. In this ta- 
ble, it is assumed that MACHOs constitute 100% of the 
halo mass. The probabilities were calculated using mag- 
nification patterns with three different side lengths (the 
three main columns on the right) for several combinations 
of MACHO mass and quasar size (indicated in the two 
columns on the left). No values are given where the pa- 
rameters were beyond the dynamical range of the simu- 
lations. The statistical uncertainties arc given in brackets 
(0.0 is given where the uncertainty was below the round- 
ing precision). Probabilities above 95% are highlighted in 
grey - the respective parameter pairs are ruled out by the 
observations at the 95% level. 



Macho Quasar 
mass size 
(M ) (cm) 



10" 



10" 



20 



pattern side length 
(Einstein radii) 

200 2000 



10 14 
3 x 10 14 

10 15 
3 x 10 



ir, 



10 14 
3 x 10 14 

10 15 
3 x 10 



ir, 



58.5(1.6) 
60.1(1.7) 
62.7(1.4) 
55.8(1.0) 

89.8(1.0) 
90.9(0.9) 
89.8(0.7) 
66.9(1.3) 



io~ 3 10 14 


99.7(0.1) 




99.8(0.0) 


3 x 10 14 


99.7(0.1) 




99.7(0.0) 


10 15 


97.1(0.6) 




96.8(0.0) 


3 x 10 15 


37.1(4.7) 




51.6(0.6) 


10~ 4 10 14 


100.0(0.0) 


100.0(0.0) 


3 x 10 14 


100.0(0.0) 


100.0(0.0) 


10 15 


89.6(1.2) 




93.8(0.2) 


3 x 10 15 






1.1(0.3) 



10~ 5 


10 14 


100.0(0.0) 100.0(0.0) 




3 x 10 14 


100.0(0.0) 100.0(0.0) 




10 15 


58.8(2.5) 60.6(0.2) 




3 x 10 15 


0.0(0.0) 0.0(0.0) 


10~ 6 


10 14 


100.0(0.0) 100.0(0.0) 




3 x 10 14 


98.0(0.4) 98.1(0.0) 




10 15 


2.4(0.3) 3.2(0.2) 




3 x 10 15 


0.0(0.0) 


io- 7 


10 14 


100.0(0.0) 




3 x 10 14 


58.4(1.3) 




10 15 


0.0(0.0) 




3 x 10 15 


0.0(0.0) 


10~ 8 


10 14 


98.4(0.2) 




3 x 10 14 


0.0(0.0) 




10 15 


0.0(0.0) 




3 x 10 15 





Table 2. Same as Table [ijfor the case where 50% of the 
halo mass is contained in MACHOs. 



Macho Quasar 
mass size 
(M Q ) (cm) 



10" 



10" 



10" 



10" 



20 



pattern side length 
(Einstein radii) 

200 2000 



10 14 
3 x 10 14 

10 15 
3 x 10 15 

10 14 
3 x 10 14 

10 15 
3 x 10 15 



53.1(2.4) 
51.9(0.9) 
53.5(1.0) 
50.1(0.8) 

85.5(0.5) 
86.1(0.6) 
86.0(0.5) 
65.1(1.7) 



10 -3 


10 14 


98.9(0.1) 


99.0(0.0) 




3 x 10 14 


99.1(0.1) 


98.9(0.0) 




10 15 
3 x 10 15 


96.5(0.4) 
39.7(9.3) 


96.0(0.0) 
36.1(0.8) 


IO" 4 


10 14 


100.0(0.0) 


100.0(0.0) 




3 x 10 14 


100.0(0.0) 


100.0(0.0) 




10 15 
3 x 10 15 


90.9(2.1) 


89.5(0.2) 
0.0(0.0) 


10~ 5 


10 14 
3 x 10 14 




100.0(0.0) 
100.0(0.0) 




10 15 
3 x 10 15 




40.2(2.6) 
0.0(0.0) 


IO -6 


10 14 
3 x 10 14 

10 15 
3 x 10 15 




100.0(0.0) 
94.4(1.0) 
0.3(0.4) 



10 14 
3 x 10 14 

10 15 
3 x 10 15 

10 14 
3 x 10 14 

10 15 
3 x 10 15 



100.0(0.0) 
100.0(0.0) 

39.4(0.2) 

0.0(0.0) 

100.0(0.0) 
93.4(0.1) 
0.3(0.0) 
0.0(0.0) 

100.0(0.0) 

25.7(2.4) 

0.0(0.0) 

0.0(0.0) 

83.1(2.6) 

0.0(0.0) 

0.0(0.0) 



ulations we find that this rules out a dominant popula- 
tion of compact objects in the halo of the lensing galaxy 
from 10 -3 Mq down to IO -5 M for quasar sizes below 
3 x 10 14 cm. These limits are almost independent of the 
fractions of the halo mass contained in compact objects. 

These results are consistent with the 2nd year results 
from the MACHO microle nsing search towards the Magel- 
lanic clouds (Alcock et al. 1997 ). The most probable mass 
of MACHOs in the Milky Way that emerges from their 
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Table 3. Same as Table [j] for the case where 25% of the 
halo mass is contained in MACHOs. 



MACHO Quasar 



pattern side length 



Q0957+561. Similar to us, they cannot find evidence for 
microlensing variations on short time scales. 

Along the same lines, our results are consistent with 



mass 
(Me) 


size 
(cm) 


(Einstein radii) 
20 200 2000 


lCr 1 


10 14 
3 x 10 14 

10 15 
3 x 10 15 


35.8(1.6) 
36.2(1.6) 
37.1(1.8) 
36.3(1.7) 




ltr 2 


10 14 
3 x 10 14 

10 15 
3 x 10 15 


74.1(1.1) 
74.7(1.2) 

54.9(2.6) 




ltr 3 


10 14 


96.2(0.5) 


95.9(0.1) 




3 x 10 14 

10 15 
3 x 10 15 


96.5(0.5) 
16.0(8.2) 


95.9(0.1) 

GO 1 (C\ 1 \ 

20.2(0.3) 


1(T 4 


10 14 


99.9(0.1) 


99.9(0.0) 




3 x 10 14 


99.8(0.1) 


99.9(0.0) 




10 15 
3 x 10 15 


73.8(7.8) 


81.5(0.2) 
0.0(0.0) 


ltr 5 


10 14 
3 x 10 14 




100.0(0.0) 100.0(0.0) 
99.9(0.0) 99.8(0.0) 




10 15 
3 x 10 15 




15.4(2.7) 18.1(0.2) 
0.0(0.0) 0.0(0.0) 


ltr 6 


10 14 
3 x 10 14 

10 
3 x 10 15 




100.0(0.0) 100.0(0.0) 

77.2(4.1) 80.7(0.2) 
0.0(0.0) 0.0(0.0) 
0.0(0.0) 


ltr 7 


10 14 
3 x 10 14 

10 15 
3 x 10 15 




100.0(0.0) 

7.2(1.2) 
0.0(0.0) 
0.0(0.0) 


1(T 8 


10 14 
3 x 10 14 

10 15 
3 x 10 15 




42.8(2.7) 

0.0(0.0) 

0.0(0.0) 



study lies around half a solar mass, which is not excluded 
by our study. In fact, such large masses are not yet probed 
by the observational data set we used (K97) because mi- 
crolensing effects by objects with masses of about one solar 
mass are only becoming observable on time-scales of sev- 
eral years. In a different data set on Q0957+561 (Schild 



1996a) there were indications for some microlensing action 
which can be produced by stars in this mass range. Very 
recently, Pelt et al. (1998) analysed a large data set on 



the results of Lewis & Irwin (|996|) who found MACHO 
masses to be in the range from 0.1M Q up to 10.0 M© from 
microlensing in the quasar Q2237+0305. 

The limits on compact objects in the halo of the lens- 
ing galaxy of Q0957+561 will, however, improve with time 
since the monitoring of Q0957+561 at the Apache Point 
Observatory (and also at other observatories) is an ongo- 
ing project. If, for example, microlensing variations above 
0.05 magnitudes are not observed for another season, the 
limits from Tables [j], || and || would become stronger and 
we could possibly rule out compact objects with masses 
as high as 0.1M© or more. 

In quasar microlensing, not only the mass of MACHOs, 
but also the size of the quasar enters the calculations. 
With the current data, however, we cannot constrain the 
quasar size very much. To do this, one would need to find 
characteristic events in the difference light curve of the two 
quasar images. If, for example, the peak in the difference 
light curve (Schild |1996b| , K97) in Fig. [j] around day 750 
is real, it would certainly be a valuable constraint for both 
the masses of MACHOs and the quasar size. 

Finally, in magnification patterns with side lengths of 
several hundred Einstein radii we find coherent structures 
in the caustic network on scales of 200 Einstein radii or 
more. In the lensing galaxy of Q0957+561 these scales cor- 
respond to physical scales of the order of l-7w jf- hr 1 ! 2 



pc 

for MACHOs of mass M. On these scales the dark objects 
might not be distributed completely randomly anymore. It 
would thus be interesting to investigate further the effect 
of clustering of MACHOs on the microlensing properties 
of galaxies. 
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